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Abstract -For a quantum many-body system, the direct population of states of double-excitation 
character is a clear indication that correlations importantly contribute to its nonequilibrium prop- 
erties. We analyze such correlation-induced transitions by propagating the nonequilibrium Green's 
functions in real-time within the second Born approximation. As crucial benchmarks, we compute 
the absorption spectrum of few electrons confined in quantum wells of different width. Our results 
include the full two-time solution of the Kadanoff-Baym equations as well as of their time-diagonal 
limit and are compared to Hartree-Fock and exact diagonalization data. 
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Introduction. — The prediction and interpreta- 
tion of electronic excitations and photoeniission (absorp- 
tion) spectra [1] is becoming more and more vital for 
jthe study of (sub) femtosecond processes in atoms [2], 
(bio)molecules [3], solids [4], and nanoscale materials [5]. 
Thereby, in many finite systems, particular importance 
must be placed on double-excitations (DEs) as they fa- 
cihtate correlated excitation pathways which considerably 
enrich the nonequilibrium properties. The final state of 
such a transition is called a doubly-excited state, i.e., it 
is one which has dominant DE character but is not neces- 
sarily well described by a single doubly-excited Slater de- 
terminant, e.g., [6]. Prominent examples include the mul- 
tiple electron-hole pair generation in semiconductors [7] 
and the autoionization of atoms which involves the inter- 
mediate population of shake-up resonance states prior to 
fragmentation [8]. 

Aside from configuration interaction (CI) and exact 
diagonalization (ED) methods^ [9], the calculation of 
electronic excitations is mainly based on time-dependent 
density functional theory (TDDFT), e.g., [1, 6], and 
many-body perturbation theory (MBPT) using Green's 
functions, e.g., [1, 10] and references therein. In both 
approaches, the central quantity is the real-time, re- 



(2),R, 



tarded two-particle response function [11,12] Xij u (^7 1') — 

^Pij i^)/^fik \^')' where p'^^ denotes the one-particle re- 
duced density matrix (IpRDM), /^^^ is an external per- 
turbation, and the poles of ^('^'•^ in frequency space refer 
to the excitation energies. The response function satis- 
fies a Bethe-Salpeter equation (BSE) [10] which contains 
a four-point integral kernel. For this reason, solutions are 
computationally challenging and often obtained at the ex- 
pense of full frequency dependence and self-consistency, 
see ref. [1] for an overview. In TDDFT, the equations 
are simpler (being of two-point-type) but depend on the 
(generally unknown) exchange-correlation functional. 

In this Letter, we apply nonequilibrium real-time 
Green's functions techniques to compute the excitation 
properties — in particular the DEs — of a correlated inho- 
mogeneous quantum system in equilibrium. As shown in 
refs. [11, 13], this allows to avoid the direct solution of 
the BSE for ;^(2),R_ Instead, we compute the first-order 
variation of the IpRDM, 



df^^(t) = -isg^^{t,t^ 



-i\im6gl-\t,t + e), (1) 



in terms of the nonequilibrium Green's function (NEGF) 
Q'^^^ [14] following a suitable external excitation^ f^^H^)- 



^CI and ED scale exponentially with system size and, hence, are ^The excitation spectrum is then obtained by inversion of the 

unfavorable for the description of large systems. linear response relation. 
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The key advantage of this method is that rather sim- 
ple conserving approximations for the propagation of the 
NEGF (such as the second Born (2B) approximation) 
translate into high-level approximations for ;^(2)'^ that 
obey the relevant sum rules [11]. This concept was ap- 
plied to excitons in optically excited semiconductors [15], 
the dynamic structure factor of the correlated electron 
gas [11] and to the absorption spectrum of small atoms 
and molecules [16,17]. 

A first extension of this real-time approach to DEs has 
been presented in ref. [12] focusing on the spectral and 
response function of small Hubbard nanoclustcrs in the 
modcrate-to-strong coupling regime. Yet no unambigu- 
ous criterion how to identify DEs in the spectrum has 
been given, and the accuracy and scope of validity of 
the used 2B for DEs has remained open. In this Letter, 
we answer these open questions by systematically analyz- 
ing electronic excitations in a few-particle quantum well 
(QW) structure as function of interaction strength (con- 
trolled by the QW width). We propagate the two-time 
NEGF fully including memory effects [18-20]. In order 
to achieve sufficiently long simulation times, we also con- 
sider the time-diagonal limit by employing the generalized 
Kadanoff-Baym ansatz (GKBA) [21]. 

Theory. — We start from a generic A^-electron Hamil- 
tonisin H^^^ ~ Ha +Hr composed of a non-interacting 
(a) and an interacting part (/3). The solution of the 
Schrodinger equation, (H^^^ — Em )|^m ) = 0, provides 
the complete spectral information including the ground 
(excited) state with energy Eq (£'m>o)- -f*-"" ^ pertur- 
bation mediated by an operator F^^' , and for the system 



initially prepared in the state [^ 



W\ 



the full excitation 



spectrum, in the linear regime, consists of discrete peaks 



at frequencies LOm 
sition moment /„ 
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for which the tran- 



IS non-zero. 



Exact diagonalization. For a given total spin 5* and 
projection Ms, the wave function |^m ') has the rep- 
resentation, 1*^^) = A^_^^{\^pin^)\S,Ms)}, where aL^' 

accounts for antisymmctrization and \tprn ) describes 
the spatial part. If the Hamiltonian H^^^^ is diagonal- 
ized in a product basis \P^ 
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tN 



^(1) 



n =1 \€') of one- 
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iN)^ 
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„(JV) 
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V- - 1 ci"',^ , \Py'i ) form eigenvectors of 



the M^ X M^ matrix iPy^', \H^^MPr'., ). This is the 
procedure of the ED technique, where typically M ^ N 
guarantees appropriate convergence. Further, as gener- 
ally not all M^ eigenvectors allow for the construction of 
a completely antisymmetric state, the relevant ground and 
excited state energies Em form an eigenvalue subset. 

Approximate excitation level (AEL). In many cases, 
the spectrum of the system can be classified according to 
the number of electrons that take part in the transition 



m — > n. Such a number is, e.g., provided by the AEL [22], 



A, 



= 7:Tr 



In im 



(2) 



Ji) 



',W\/„/,Wi 



with3 p\^> ^ Tr2...jv|V'm OiV"™ Y, where both IpRDMs, 
Pm and pii , are expressed in the natural orbital basis 
that diagonalizes the initial state density matrix. While 
transitions with an AEL close to unity are referred to as 
single-excitations (SEs), an A„_j.„ « 2 (3 or 4 etc.) in- 
dicates states of essential double- (triple- or quadruple- 
etc.) excitation character. We note, that an integer 
AEL is only obtained in an effective single-particle ap- 
proach using, e.g., the occupied and virtual Kohn-Sham 
or Hartrec-Fock (HF) orbitals |Q ') of the ground state. 
In this frozen-orbital (FO) or Koopmans' approxima- 
tion [9], one or more electrons arc promoted into vir- 
tual orbitals forming single singly- or multiply-excited 
Slater determinants. The corresponding ground (excited) 
state energy is approximate due to the neglect of or- 
bital relaxations and correlations^: £",, — "^' ^ 



y-' t(2) 



Si7 ^i" • where M; denote the orbital en- 



ergies and jj;f = {ii\hf^\jj) {K'>f = {ij\'hf^\ji)) are the 
Coulomb (exchange) integrals with {ij\hp '\kl) = h\ ^ ^.j = 
//d^r dV C,(rl Q{r') hf\f- r') a{r) 0(r '). 

Nonequilihrium Green's functions. Since CI and ED 
can only handle small systems and pure states, in the fol- 
lowing, we resort to a NEGF approach. As discussed in 
the introduction, eq. (1) involves the NEGF g2^\t,t') = 



— ^ {Tc Ci{t) c'Jf)) of the perturbed system, where c] (ci) 
denote fermionic creation (annihilation) operators, i in- 
volves spin and orbital degrees of freedom and ^| (t, t') = 

Ocit - t')gi^^->it,t') - ec{t' - 1) g\f''^{t,t'). The negf 

is defined on the round-trip Keldysh contour C [23] with 
the contour-ordering operator Tq and obeys the two-time 
(non-Markovian) Kadanoff-Baym equations (KBEs) [14], 



d 



.(1) 



(1)^ 



{i^n7'5.fe-<kW}C(*'^') 



(3) 



(1)^ 



zWi 



6c{t-t')S,,+ di^l>(t,i)g'^>{i,t') 



particle states \(f>\ ), the coefficients in the expansion , (i) 



together with the adjoint equation with t ^ t' . Here, 



(t) = {i\ha {t)\j) denotes the one-particle energy and 
summation over k is implied. A many-body approxima- 
tion (MBA) for the two-time self-energy E^^' is obtained 
by MBPT as functional of the NEGF and the interac- 

(2) 

iiki- ^^ ^^ approximation, one obtains 

(1),HF/ 



tion potential h 

.(1),HF 
'^3 



ny^- {t,t') ^6cit- 1') i:Yf^' (t) with. 



.(1),HF 



y 



(t) ^-ih{h' 



(2) 
I3,ijkl 



f^ikA5l^(t^t) 



(4) 



■'The prime indicates a second set of coordinates. 
■'Here, the prime impHes summation over only those orbitals that 
contribute in the determinant. 
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whereas in the 2B approximation one adds, 



.(1),2B 



(t,t') 



,(2) 



-(ih) 'l^^ifems {"/3^;j>n 



.(2) 






7(1)^ 



xgi\>{tX)GmtX)G^\\t',t) 



W( 



(5) 



Both expressions conserve particle number, total momen- 
tum and energy. Computationally, the time-non-local, 
second-order self-energy (5) involves a large number of 
0{n^) operations at basis dimension n^,. However, its 
evaluation can be drastically simplified by using the finite 
element-discrete variable representation^ (FEDVR) that 
leads to a scaling behavior of 0{n^), see ref. [19]. 

In this Letter, we apply two variants of solving the KBEs 
within approximations (4) plus (5) for Y.'-^K 

I. We, without any further approximations, propagate 
the NEGF in the double-time plane under the pres- 
ence of the full memory kernel (r.h.s. of the KBE). 
This gives rise to large memory requirements that 
limits the propagation time, though efficient code par- 
allelization and FEDVR-type representations consid- 
erably extend the range of applicability, see ref. [20] . 
To overcome this limitation, 

II. we employ the GKBA^ [21], g(i)^^(t,i') « 

±{gW'^{t,t') p(^^^^{f) - p^^'>^^{t)g(^'>^^{t,t')}, with 



:(i) = ±ifig(i)'^(t,i), where the Green's func- 



tions are reconstructed from their values on the time 
diagonal. This procedure is known to yield reliable re- 
sults for weak to moderate coupling, e.g., [15,25]. In 
the following, the retarded (advanced) function tj(i)'^ 
(g(i).A-j jg taken at the HF level allowing to signifi- 
cantly reduce the computational effort in comparison 
to full 2B (case I) and allowing for longer propagation 
times combined with larger basis dimensions. 

For details on the propagation of the NEGF in either cases, 
the reader is referred to refs. [18,20,26]. 

To determine the excitation spectrum, we start the 
propagation at a time Iq, when the system is in the ground 
state^. In the early stage of the evolution, the many-body 
system is perturbed in the form h^^j + fofij^{t — to). 
Complying with the constraint of sufRciently small ampli- 
tudes /o, we then compute the linear response Sp^, (t) = 



-i{gl'^'^{t,t) 



Gij {to, to)}, and the Fourier transform 



of (/'^^)(t) = /oTr[(5p(i)(t)/(i)] yields the absorption spec- 
trum and the excited state energies. 

Model. — As test system, we consider A^ = 4 elec- 
trons (charge e, effective mass m* , positions Vi) with sin- 
glet spin configuration \S,Ms) — |0,0) in a quantum 
well (QW) potential of width L. Neglecting the lateral 



^Here, hg ^ ^.j becomes diagonal in, both, ij and kl. 

®The GKBA retains the conservation properties [24]. 

^While, in case I, the initial state is obtained from the solution 
of the Dyson equation, see [19], for the GKBA (case II), we use 
adiabatic switching [27] for its generation. 



Table 1: Four-electron QW with \S, Ms) = [0, 0) (singlet state) 
at coupling A* = 5: Ground-state excitation energies cuo^n ob- 
tained from exact diagonalization (ED) and the frozen-orbital 
(FO) approximation. Ao^n denotes the AEL, eq. (2). DEs of 
dipole character are underlined. See also fig. 1 (lower panel). 



n 


^Jf5n {Eo) 


Al^n 


Character 


OjE^n {E5) 


1 


1.1445 


1.079 


dipole 


1.0616 


2 


1.5478 


1.088 


non-dipole 


1.5218 


3 


2.1118 


1.842 


non-dipole 


2.1491 


4 


2.5559 


1.150 


non-dipole 


2.4841 


5 


2.6769 


1.870 


dipole 


2.5841 


6 


2.9794 


1.145 


dipole 


2.9499 


7 


3.0394 


1.969 


non-dipole 


3.0646 


8 


3.6314 


1.944 


dipole 


3.5629 


9 


3.9928 


1.997 


non-dipole 


3.8556 


10 


4.2129 


1.956 


non-dipole 


4.0034 


11 


4.3150 


1.119 


dipole 


4.2735 


12 


4.5544 


2.048 


dipole 


4.4921 


13 


4.7619 


1.122 


non-dipole 


4.7512 


14 


4.9300 


1.998 


non-dipole 


4.9824 


15 


4.9575 


2.830 


dipole 


4.9350 


16 


5.4099 


2.797 


non-dipole 


5.3985 


17 


5.4250 


2.034 


non-dipole 


5.3527 


18 


5.5413 


2.019 


dipole 


5.4311 


19 


5.8032 


2.038 


dipole 


5.6819 


20 


5.8882 


2.153 


non-dipole 


5.7953 


21 


6.0119 


1.981 


dipole 


5.9254 


22 


6.3480 


2.055 


non-dipole 


6.2934 


23 


6.4185 


2.839 


non-dipole 


6.3540 


24 


6.4711 


1.096 


non-dipole 


6.4505 


25 


6.7825 


2.802 


non-dipole 


6.7618 


26 


6.8345 


2.018 


dipole 


6.7711 


27 


6.8966 


2.941 


dipole 


6.8230 


28 


6.9333 


1.104 


dipole 


6.9334 


29 


7.1216 


2.108 


non-dipole 


7.0143 


30 


7.2410 


2.912 


dipole 


7.2248 


31 


7.4516 


2.035 


non-dipole 


7.2317 


32 


7.5683 


1.986 


dipole 


7.5295 


33 


7.7324 


3.876 


non-dipole 


7.7465 


34 


7.7932 


2.150 


dipole 


7.7255 


35 


7.9891 


2.031 


non-dipole 


7.8790 



electron motion, the Hamiltonian in units of the confine- 
ment energy^ Eq ~ h? / [m* L"^) reads in vertical direc- 
tion 7^W(t) = E^i h!-^\z,,t) + y Ef<.<, hf{z, - z,), 

-hf'Kz) Sit-to) 



withri = (0, 0, ZiL), ha {z,t) — 

:(2) 



1_9_ 
'2 dz^ 



and hy{z- z') = [{z - z'f + k2]-i/2_ The coupling 
strength A* = Lja}^ = e^m*Ll{A.'KeoS*f?) with effective 
Bohr radius Oq and (material) dielectric constant (e*) £o 
defines the relative interaction strength between the elec- 
trons. While A* — > represents the ideal quantum regime, 
the limit A* ^ oo leads to quasi-classical, Wigner-crystal 
behavior. In a GaAs heterostructure (m* = O.OGTwe, 

*i?g7r^/2 is the ground-state energy of a single-electron QW. 
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Fig. 1: (Color online) Ground-state excitation spectrum of the four-electron QW at A* = 5. Lower panel: approximate excitation 
level, eq. (2). The red open (closed) dots indicate states that are accessible from the ground state by (non-)dipole transitions. 
The black squares show the excitation energies in FO approximation (integer AEL). All dipole-allowed SEs and DEs (cf. the 
black vertical lines) are numbered according to table 1. Upper panel: Dipole excitation spectrum \f^^' {uj)\^ as obtained from 
the time-dependent Schrodinger equation (TDSE) and from the Kadanoff-Baym equations in HF (black dashed line) and 2B 
approximation (blue dots). The labels a-e mark prominent deviations of the 2B result from HF and TDSE, respectively. 



e* = 12.9), the moderate-to-strong coupling cases A* = 1 
(5) correspond to well widths of L = 10.2 nm (50.9 nm) 
which are realistic QW dimensions at confinement energies 
E^ = 10.956 meV (0.438 meV). Further, the soft-Coulomb 
form of h\ with regularization k, (we set k = 0.2) pre- 
vents the divergence of the two-electron integrals which are 
constructed from one-dimensional FEDVR spin-orbitals. 
In the following, we consider /^^^ as the dipole opera- 
tor, i.e., f^^'{z) = z such that /o is measured in units of 
eLSQ/EQ, where £q is the electric field strength. The cor- 
responding transition dipole moment (TDM) is defined by 



/. 



(N) 



e{ipn 



Wi 



YJi'=iZtH' 



W\ 



For transitions of dipole 



(non-dipole) type, it is |/, 



(N) 



> (I/, 



{N) 



0). 



Results. — First, we concentrate on the moderate- 
to-strong coupling case A* = 5 which is typical for sub- 
100 nm semiconductor heterostructures. In contrast to the 
ideal, non-interacting system^, where dipole selection rules 
prohibit DEs, we, here, expect doubly-excited states of 
relevant TDM. Table 1 shows all ground-state excitation 
energies a;™„ below 8E'^ as obtained from ED itemized 
by AEL and general dipole character. We observe, that 
dipole and non-dipole transitions alternate and that the 
AEL allows for a unique classification of each excited state. 
Overall, its deviation from an integer is less than 0.21. 

According to table 1, the dipole excitation spectrum 
starts with a SE (n = 1) of energy ujq^i = 1.14 _Eg fol- 
lowed by a DE (n = 5) with wfi^g = 2.68 E'jJ, a SE 
(n = 6) with w™6 = 2.80 E^, and a DE (n = 8) with 
a;™g = 3.63 E'g, etc.. The first dipole-allowed triple- 
excitation has quantum number n — 15. Further, we note, 
that the energetically lowest doubly-(quadruply-)excited 
state with n = 3 (33) is not of dipole-type, cf. also fig. 1 

®The ideal ground-state energy is ■7t^/5Eq. 



(lower panel). Moreover, table 1 includes the result of 
the FO approximation as carried out using the ground- 
state HF orbitals, see w™„. Thereby, the energies well 
approach the exact ones although FO can lead to wrong 
ordering, cf. n = \A vs. 15 (16 vs. 17), and can produce 
shifts as large as 0.2 i?Q, cf. n = 10 and 31. 

The situation for A* = 5 is visualized in the lower and 
upper panel of fig. 1 including exact and the approxi- 
mate NEGF results. From the linear-response dipole spec- 
trum as obtained from the solution of the time-dependent 
Schrodinger equation (TDSE), see the red solid line for 
|/^i^(aj)p (upper panel), we find that DEs appear indeed 
with significant TDM, cf., e.g., the peaks with n = 5 and 
8. From the solution of the KBEs (3) according to case I, 
we obtain the time-dependent HF and 2B result, which, 
respectively, are indicated by the black dashes and the 
blue dots. First, we observe that HF has only the ability 
to describe SEs, cf. the peaks for states with n = 1, 6, 
11 and 28. Thereby, the corresponding energies are more 
accurate than those in FO approximation due to the in- 
clusion of orbital relaxations. The failure of reproducing 
peaks of multiply-excited states, however, leads overall to 
a very simple spectrum. In contrast, the 2B approxima- 
tion, which accounts for more than 70% of the correlation 
energy of the QW's electronic ground state^", shows much 
more structure in the spectrum. In particular, it yields 
additional peaks that are located at energies for which we 
expect DEs, cf. the states indicated n = 5, 8 and 12. 

It is tempting to identify the additional peaks in 2B as 
the correct DEs. However, great care is needed — especially 
in the context of additional deficiencies in the 2B result 
which require explanations: (i) the relatively low excita- 



^"At A* = 5, the exact ground-state energy is 5.5278 -Bq. 
The approximate treatment yields in HF: 5.5384 Bg, in 2B 
(case I): 5.5306 E^, and in GKBA (case II): 5.5298 E^. 
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Fig. 2: (Color online) Ground-state dipole excitation spectrum 
of tiie four-electron QW as function of the coupling strength 
A*. HF: green dashed lines, GKBA (case II): gray-scale density 
plot, full 2B (case I): blue dots, ED: red solid lines, and FO for 
DEs with j4o->n = 2: black squares. Note, that, for ED, only 
SEs and DEs are shown. While labels a-c mark the avoided 
crossings in the GKBA result, labels d-f indicate excitations 
which are not attributable to DEs. 



tion strength for the peaks at n = 5 and 8 relative to 
those at n = 6 and 11 (compare with the TDSE data), (ii) 
the observation of energy shifts and splittings for singly- 
excited states which are well covered even by the HF result 
(see labels a, b), and (iii) the presence of peaks at low and 
high energies which cannot be attributed to DEs of the 
system (cf. labels c, d and e). 

With the intention to explain these shortcomings, we 
systematically investigate the excitation spectrum as func- 
tion of A*. Focusing on SEs and DEs, fig. 2 shows how the 
spectrum changes when going from the non-interacting 
system to the moderate-to-strong coupling case A* = 6. 
In HF approximation, the SEs with n = 1, 6, 11 and 28 
are very well described over the whole A*-range, see the 
green dashed lines which practically lie on top of the ex- 
act result (red lines). For small A* < 1, the same is found 
for the GKBA result (gray-scale density) such that, in the 
limit A* -> 0, for, both, HF and GKBA, the SE ener- 
gies exponentially converge towards the ones of the ideal 




Fig. 3: (Color online) Four-electron QW as in figs. 1 and 2. 
Shown is the absolute square of the TDMs for the five energet- 
ically lowest excited states as function of the coupling strength 
A*; states with n = 1, 6, 11: SEs, n = 5, 8: DEs. GKBA: black 
lines, 2B: blue dots, and ED: red solid lines. For A* > 2, 
the TDM for the energetically lowest DE (label a) indicates a 
transition from a DE to a SE, compare with fig. 2. 



system. For doubly-excited states, we find a completely 
different behavior. Comparing GKBA to ED and to FO 
(black squares), we observe that the DE energies as de- 
scribed by GKBA, vary too strongly with the coupling 
strength though the correct A* — 0-limit is reached. Even 
"avoided" crossings of the A^-particle energy levels occur 
for A* ~ 1.75 and thereafter for A* ~ 3.5 and 5, cf. labels 
a, b and c. Such a behavior is not supported by the exact 
description^^. However, it is not a deficiency of the GKBA 
because full 2B calculations (see the blue dots which mark 
resolvable peak positions) confirm this behavior. 

To prove that the non-singly excited states found by 
GKBA and 2B are really of DE character, we, in fig. 3, 
analyze the TDMs of the five lowest excited states as func- 
tion of A*. First, for the three SEs below A* ~ 1, the ex- 
act TDMs are nearly constant, see the red lines indicated 
n = 1, 6 and 11. Second, for the two DEs (n = 5 and 8), 
l/oVnP is several orders of magnitude smaller for low cou- 
pling strength and increases with (A*)^. Both these prop- 
erties can be understood by standard perturbation theory 
for the four-electron QW and are well reproduced by the 
GKBA calculations, see the black lines. Beyond a coupling 
strength of A* w 0.5, however, the approximate and exact 
TDMs start to deviate considerably. Here, GKBA shows 
a complex A* -dependence whereby the energetically low- 
est (isolated) excited state with quantum number n — \ 
is less affected. We attribute this to the "avoided" energy 
levels crossings, where the energetically lowest DE takes 
over the role of the second SE (cf. curve a and compare 
with fig. 2). The same occurs for the second DE at A* ~ 4. 



'At least in the considered regime of coupling strengths. 
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Discussion. — Following figs. 2 and 3, we can clearly 
distinguish between states of different excitation character 
in the GKBA (2B) spectrum and conclude that reasonable 
results are obtained for small coupling. Moreover, we are 
now in a position to give reasons for some of the defi- 
ciencies formulated for the A* — 5-case. Obviously, the 
presence of an "avoided" crossing for the third (fourth) 
SE in the case of 2B and GKBA as function of A* explains 
the energy shift (splitting) in the corresponding excita- 
tion energy, cf. label a (b) in fig. 1. The fact that the 
2B approximation underestimates the TDMs for several 
excited states is of similar reason, see, e.g., the GKBA 
result for the regime A* > 2 in fig. 3. A further main 
results is that, due to the strong coupling dependence of 
the DE energies in 2B, the peak at w « 3.6 E'q {n — 8) 
in fig. 1 (blue dots) does not originate from the second 
but from the first DE, compare with fig. 2. In this re- 
gard, we also note that practically all energy differences 
for doubly-excited states do not exponen- 



W( 



2B 



to, 



ED 

0->n 



tially become small in the limit A* — > 0, although this is 
the case for all SEs. Instead, it seems that DEs have in 2B 
an intrinsic coupling dependence determined by perturba- 
tion character of the MBA. Furthermore, we mention, that 
the additional peaks emerging in the 2B result disappear 
for small A* but cannot be attributed to states of the real 
system^^, cf. labels c-e in fig. 1 (upper panel) and labels d 
and e in fig. 2. The study of their origin requires further 
investigations^^ — however, we report that their presence 
does not depend on the approximation level of the ground 
state a.t t — to, i.e., whether it is of HF-type or correlated. 

In conclusion, the four-electron QW is sufficiently sim- 
ple to allow for ED results. This has given us the pos- 
sibility to quantify the accuracy of the 2B approxima- 
tion with respect to important correlation features such 
as electronic DEs, which gain importance for moderate to 
strong coupling (cf. fig. 3). The spectrum follows from 
direct time-propagation of the NEGF, where the GKBA 
has allowed for a detailed analysis for a broad range of in- 
teraction strengths (tuned by the QW width). DEs have 
been identified by tracing the system behavior to the zero- 
coupling limit analyzing the functional dependence of their 
oscillator strengths. In contrast to ED, we, in the 2B re- 
sult (with and without the GKBA), have found a strong 
coupling dependence of the DE energies and, as a conse- 
quence, mixings (or hybridizations) of SEs and DEs which 
show up as "avoided" crossings, cf. fig. 2. As this occurs 
in the context of well-described SEs in a HE treatment, 
great care is needed when interpreting correlated excita- 
tion pathways in MBAs beyond HE. It will be interesting 
to see, whether and how much the present behavior can 
be "repaired" by more advanced many-body self-energies 
such as the T-matrix or GW approximation. 



-'•^Notc, that these artifacts ean mimic relevant excited states, cf., 
e.g., label f in fig. 2 (A* = 5) which refers to the first DE. 

^^Similar effects have been observed in the 2B treatment of exci- 
tations in small Hubbard chains, cf. ref. [26]. 
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